Can the nuclear symmetry potential at supra-saturation densities be negative? 
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In the framework of an Isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport 
model, for the central ^''^Au+^^^Au reaction at an incident beam energy of 400 MeV/nucleon, 
effect of nuclear symmetry potential at supra-saturation densities on the pre-equilibrium clusters 
emission is studied. It is found that for the positive symmetry potential at supra-saturation den- 
sities the neutron to proton ratio of lighter clusters with mass number A < 3 {{n/p)A<3) is larger 
than that of the weighter clusters with mass number ^ > 3 {{n/p)A>3), whereas for the negative 
symmetry potential at supra-saturation densities the {n/p)A<3 is smaller than that of the {n/p)A>3- 
This may be considered as a probe of the negative symmetry potential at supra-saturation densities. 
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Recently the studies of the density-dependent nuclear 
symmetry energy, which is crucial for understanding 
many interesting issues in both nuclear physics and as- 
trophysics [i|-@|> have attracted much attention |9l-lll|. 
The high density behavior of the symmetry energy, how- 
ever, has been regarded as the most uncertain property 
of dense neutron- rich nuclear matter |12l . Il3| . Many mi- 
croscopic and/or phenomenological many-body theories 
using various interactions predict that the symmetry en- 
ergy increases continuously at all densities. On the other 
hand, other models predict that the symmetry energy 
first increases to a maximum and then may start de- 
creasing at certain supra-saturation densities. Thus, cur- 
rently the theoretical predictions on the symmetry energy 
at supra-saturation densities are extremely diverse. To 
make further progress in determining the symmetry en- 
ergy at supra-saturation densities, what are most crit- 
ically needed is some guidance from dialogues between 
experiments and the predictions of heavy-ion collisions 
transport models, which have been done extensively in 
the studies of nuclear symmetry energy at low densities 

While studying the symmetry energy by using heavy- 
ion collisions, the related input in an IBUU transport 
model is actually the symmetry potential, which is a more 
complete information than the density dependence of the 
symmetry energy at zero temperature calculated with 
the same (mean-field) approximation. Unfortunately, the 
symmetry potential is also rather uncertain, which can be 
positive or negative at supra-saturation densities [3, US~ 
|24| . To study the symmetry energy at supra-saturation 
densities in the framework of an IBUU model, one has to 
firstly determine the symmetry potential at higher den- 
sities. The symmetry potential, its positive or negative, 
thus urgently needs to be solved. In this paper, we use 
the asymmetry of cluster emission in dynamical simula- 
tion to study the symmetry potential at supra-saturation 
densities. 

The non-equal partition of the system's isospin asym- 
metry with the gas phase being more neutron-rich than 
the liquid phase has been found as a general phenomenon 
using essentially all thermodynamical models and in 
simulations of heavy-ion reactions [2, l25l - [28| . But all 



these studies are for low energy density nuclear matter, 
in which the sub-saturation symmetry energy/potential 
dominates the isospin fractionation. For supra-saturation 
symmetry energy/potential's studies, one needs to study 
the isospin fractionation of higher energy density nu- 
clear matter. Such matter is explored in the first step 
of high energy heavy-ion collisions, where time depen- 
dence and out-of-equilibrium effects play a dominant 
role. Therefore, to study the high density behavior of 
nuclear symmetry energy/potential, one needs to base 
on the transport model, to study particle emission by 
using relative high incident beam energy of heavy-ion 
collisions and compare with the experimental data. In 
the present study, using an isospin and momentum- 
dependent transport model IBUU, as an example, we 
studied the neutron-rich reaction of ^^^Au-t-^^^Au at 
a beam energy of 400 MeV/nucleon with the positive 
and negative symmetry potentials at supra-saturation 
densities while keeping the low density symmetry en- 
ergy/potential fixed. We compute the average isospin 
ratio of clusters of size larger (smaller) than A=4 (corre- 
sponding mass number of a particle), which will be called 
{n/p)A>3 {{n/p)A<3) in the following. We find that for 
the positive symmetry potential at supra-saturation den- 
sities the {n/p)A<3 is larger than that of the {n/p)A>3, 
whereas for the negative symmetry potential at supra- 
saturation densities the {n/p)A<3 is smaller than that of 
the {n/p)A>3- 

The isospin and momentum-dependent mean field po- 
tential (MDI) used in the present work is [29] 
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where S — (p„ — Pp)/p is the isospin asymmetry of the 
nuclear medium. In the above t = 1/2 (—1/2) for neu- 
trons (protons) and t j^ t'; a = 4/3; /t(i", p) is the phase 



space distribution function at coordinate r and momen- 
tum p. The parameters A„(a:), A;(x), B, Ct-^t-, Ct-^t-/ and 
A were obtained by fitting the momentum-dependence 
of the U{p,S,p,T,x) to that predicted by the Gogny 
Hartree-Fock and/or the Brueckner-Hartree-Fock (BHF) 
calculations [30|, the saturation properties of symmet- 
ric nuclear matter and the symmetry energy of about 30 
MeV at normal nuclear matter density po = 0.16 fm^^ 
|19l |29| . The incompressibility Kq of symmetric nuclear 
matter at po is set to be 211 MeV consistent with the 
latest conclusion from studying giant resonances [3ll - t33| . 
The parameters Au{x) and Ai{x) depend on the x pa- 
rameter according to 
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The variable x is introduced to mimic different forms 
of the symmetry energy/potential predicted by various 
many-body theories without changing any property of 
the symmetric nuclear matter and the symmetry energy 
at normal density pQ. The last two terms in Eq. ([T]) 
contain the momentum-dependence of the single-particle 
potential. The momentum dependence of the symmetry 
potential stems from the different interaction strength 
parameters C-r,T' and Cr.r for a nucleon of isospin r inter- 
acting, respectively, with unlike and like nucleons in the 
background fields. More specifically, we use CunUke = 
-103.4 MeV and Cuke = -11.7 MeV. With these pa- 
rameters, the nucleon isoscalar potential estimated from 
t^isoscaiar ~ (t^n+t^p)/2 agrecs with the prediction of vari- 
ational many-body calculations for symmetric nuclear 
matter [T9I, ^29, Isi] , the BHF approach 21,'3^,'35!l includ- 
ing three-body forces and the Dirac-Brucckner-Hartree- 
Fock (DBHF) calculations (Sg in broad ranges of density 
and momentum. And the corresponding pressure of sym- 
metric matter is consistent with the experimental limits 
[37| . The corresponding isovector (symmetry) potential 
can be estimated from C/sym ~ {Un — Up)/25. With differ- 
ent X parameters at normal nuclear matter density po ■, the 
symmetry potential, as shown in Fig. [U agrees very well 
with the Lane potential extracted from nucleon-nucleus 
and (n,p) charge exchange reactions available for nucleon 
kinetic energies up to about 100 MeV [38]. At supra- 
saturation densities we can see that the stiff symmetry 
energy (x — —2, shown in Fig. [5]) corresponds to posi- 
tive symmetry potential while the soft symmetry energy 
corresponds to negative symmetry potential. 

According to essentially all microscopic model calcula- 
tions, see e.g., [30, Is^, the EOS for isospin asymmetric 
nuclear matter can be expressed as 



E{p, 5) = E{p, 0) + E,y^{p)S^ + OiS""), 
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where E{p, 0) is the energy per nucleon of symmetric nu- 
clear matter, and i?sym(p) is the nuclear symmetry en- 
ergy. With the single particle potential Eq.(IT]), the sym- 
metry energy can be written as [39| 
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FIG. 1: The symmetry potential as a function of energy for 
different density values in the MDI interaction with x — 1 
and X = —2. The experimental data [3g are also shown. 
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where pf = ?i(37r^|)^'"^ is the Fermi momentum for sym- 
metric nuclear matter at density p. In the mean-field ap- 
proximation, the symmetry energy Esym {p, x) is actually 
decided by its potential part EPyj^{p,x), Esymip,x) = 
g^(3^2^/2)2/3 + EPyJp,x) and the EPy^{p,x) is di- 
rectly related to the symmetry potential (C/n — Up)/ 26 
[19|. Shown in Fig. [2] is the density-dependent symmetry 
energy. At sub-saturation densities, recent studies have 
constrained the symmetry energy around 31.6{p/ pq)^'^^ 
(corresponding a; = 0) [3, [l^ |40, |4i|. So in the present 
work, we use the symmetry energy/potential correspond- 
ing a; = as our choice at sub-saturation densities. 
At supra-saturation densities, because the symmetry en- 
ergy/potential is very uncertain, we use the positive 
(x — —2) and negative {x = 1) symmetry potential as 
two extreme cases (corresponding super-stiff and super- 
soft symmetry energy at supra-saturation densities, re- 
spectively). 

In the IBUU transport model, the initial neutron and 
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FIG. 2: The symmetry energy as a function of density in the 
MDI interaction with different x parameters. The symmetry 
energy at sub-saturation densities was roughly constrained 
from recent studies 0, E, 113, Ii3l • 




20 40 60 80 100120140160 20 40 60 80 100120140160 



E (MeV/nucleon) 



FIG. 3: Neutron to proton ratio n/p as a function of nucleonic 
kinetic energy of clusters with A < 3, A > 3 from the central 
^^^Au+^^^Au reaction at a beam energy of 400 MeV/nucleon 
with the positive {x — —2) and negative {x — 1) symmetry 
potentials at supra-saturation densities but keeping the sym- 
metry potential at sub-saturation densities fixed (a; = 0). The 
number 1.49 denotes the reaction system's n/p. 



proton density distributions of the projectile and target 
are obtained by using the Skyrme-Hartree-Fock theory. 
The isospin-dependent in-medium nucleon-nucleon (NN) 
elastic cross sections from the scaling model according to 
nucleon effective masses are used. For the inelastic cross 
sections we use the experimental data from free space NN 
collisions since the in-medium inelastic NN cross sections 
are still very much controversial. The total and differen- 
tial cross sections for all other particles are taken either 
from experimental data or obtained by using the detailed 
balance formula. The isospin dependent phase-space dis- 
tribution functions of the particles involved are solved by 
using the test-particle method numerically. The isospin- 
dependence of Pauli blockings for fermions is also con- 
sidered, for more details we refer the reader to Ref. [10| . 
The nuclear liquid-gas phase transition in dilute asym- 
metric nuclear matter is studied recently [42l |. It argued 
that the neutron to proton ratio of the gas phase be- 
comes smaller than that of the liquid phase for energetic 
nucleons and the gas phase is still overall more neutron- 
rich than that of liquid phase. Noticing that the re- 
cent comparisons between IBUU calculations and FOPI 
data favor a rather soft symmetry energy (a; = 1) (corre- 
sponding a negative potential at supra-saturation densi- 
ties) at supra-saturation densities [9'] and the progress of 
the studies of nuclear symmetry energy at sub-saturation 
densities [3, [la |40, |4l[ , in this work we studied the n/p 
of clusters emitted in the central ^^^Au-|-^^^Au reaction 
at a beam energy of 400 MeV/nucleon with the posi- 
tive {x — —2) and negative {x = 1) symmetry poten- 
tials at supra-saturation densities but keeping the sym- 
metry potential at sub-saturation densities fixed (a; = 0). 
The maximal compressed density reached in this reac- 
tion is about 2.5po [9|. Shown in Fig. |3]is n/p of clus- 



ters with A < 3, A > 3 'AS a, function of nucleonic ki- 
netic energy. For the negative symmetry potential, neu- 
trons trend to being attracted by the symmetry potential 
and protons trend to being repelled during isospin frac- 
tionation. Thus n/p of weighter clusters (A > 3) are 
larger than that of the lighter clusters {A < 3). Whereas 
for positive symmetry potential, neutrons trend to being 
repelled by the symmetry potential and protons trend 
to being attracted during isospin fractionation. There- 
fore n/p of weighter clusters (A > 3) are smaller than 
that of the lighter clusters {A < 3). From Fig. |3l we 
can also see that n/p of weighter clusters {A > 3) with 
negative symmetry potential is overall larger than re- 
action system's n/p while with the positive symmetry 
potential n/p of weighter clusters {A > 3) is expected 
overall smaller than reaction system's n/p. In Fig. |31 
(n/p)A<3 and {n/p)A>3 are the neutron to proton ra- 
tios of nucleons with local densities smaller {A < 3) and 
larger {A > 3) than pc = l/8po, respectively. Chang- 
ing pc from l/5po to l/lOpo, our quantitative results 
only shift about 3%. The integral {n / p) a<3 I {n / p) a>3 
with nucleonic kinetic energy E kin I nucleon > 20 MeV 
(nucleons with kinetic energy E^m / nucleon < 20 MeV 
are mainly from cluster decays [lg|) is 1.44/1.65 for the 
negative symmetry potential and 1.48/1.35 for the pos- 
itive symmetry potential. In the standard implementa- 
tion of cluster recognition after the BUU dynamics [43] , 
a physical fragment is formed as a cluster of nucleons 
with relative momenta smaller than Pq = 263 MeV/c 
(Fermi-Momentum of normal nuclear matter) and rela- 
tive distances smaller than i?o = 3 fm (deduced by the 
uncertainty relationship of quantum mechanics). The in- 
tegral {n / p) A<z I {n / p) A>3 with nucleonic kinetic energy 



TABLE I: The clusterisation method and the associated pa- 
rameters. 



set 


n/p 


Pc = l/8po 


Po = 263 MeV/c 


Negative potential 
(x=l) 


A<3 
A>3 


1.44 
1.65 


1.44 
1.55 


Positive potential 

(x= -2) 


A<3 
A>3 


1.48 
1.35 


1.49 
1.42 



Ekin/nudeon > 20 MeV is 1.44/1.55 for the negative 
symmetry potential and 1.49/1.42 for the positive sym- 
metry potential and our results are not sensitive to the 
phase-space coalescence parameter settings (when chang- 
ing Po from 263 MeV/c to P£, » 258 MeV/c, the results 
shift not more than 0.3%). In order to show more clearly 
the above results, We present Table ID It is seen that the 
clusterisation method based on phase-space coalescence 
and the method based on a density cut-off give the same 
(n/p)A<3, but different {n/p)A>3- Specifically, the sym- 
metry energy dependence of the isospin content of large 
clusters is less pronounced with the coalescence method. 
This can be qualitatively understood. Indeed the nu- 
clear Fermi momentum of asymmetric nuclear matter is 
p^ = (1 - Sy^^po, P^ = (1 + Sy^^po. Where p^p^p^ are 
proton and neutron's Fermi-Momenta, respectively j35| . 
At densities close to saturation neutrons will then have a 
momentum that can exceed the coalescence parameter Pq 
which corresponds to symmetric nuclear matter at sat- 
uration. Therefore in average the coalescence criterium 
will be harder to fulfill for neutrons than for protons at 
the same density. This can explain why, for the nega- 
tive symmetry potential, {n/p)A>3 is lower with the co- 
alescence method than with the density cut-off method. 
The opposite effect is observed with the positive sym- 
metry potential, which leads to more neutron-rich dense 
matter. From the above we can see that although the 
two methods (the clusterisation method based on phase- 
space coalescence and the method based on a density 
cut-off) give different quantitative results, they give the 



same qualitative results, i.e., {n/p)A<3 is larger/smaller 
than {n/p)A>3 for the positive/negative symmetry po- 
tential at supra-saturation densities. 

The overall {n/p)A<3 is smaller than {n/p)A>3 with 
nucleonic kinetic energy Ekin/nucleon > 20 MeV for 
the negative symmetry potential at supra-saturation den- 
sities after isospin fractionation is different from the 
knowledge in the literatures [2, 25 28, 42] for lower en- 
ergy density nuclear isospin fractionation where mainly 
the symmetry energy/potential at sub-saturation den- 
sities works. In case experimentally in the central 
^^^Au-j-^^^Au reaction at an incident beam energy of 400 
MeV/nucleon {n/p)A<3 is overall smaller than {n/p)A>3 
with nucleonic kinetic energy Ekin/nudeon > 20 MeV, 
then a negative symmetry potential at supra-saturation 
densities is obtained. Otherwise the negative symmetry 
potential at supra-saturation densities is ruled out. This 
test in fact does not dependent on our quantitative results 
of {n/p)A<3 or {n/p)A>3- This can be a probe of the neg- 
ative symmetry potential at supra-saturation densities. 
It should be mentioned that different isospin dependent 
in-medium nucleon-nucleon cross sections may affect our 
present results [ij, E^ , such studies are planned. 

In summary, based on an isospin dependent transport 
model IBUU, n/p of the clusters as a test of negative 
symmetry potential at supra-saturation densities is stud- 
ied. We found that for the positive symmetry potential 
at supra-saturation densities the neutron to proton ra- 
tio of lighter clusters with mass number A < 3 is larger 
than that of the weighter clusters with mass number 
A > 3, but for the negative symmetry potential at supra- 
saturation densities the {'n/p)A<3 is smaller than that of 
the {n/p)A>3- This may be considered as a probe of the 
negative symmetry potential at supra-saturation densi- 
ties. 
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